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Inspired by the abundant experimental observation of axial-vector states, we study whether the observed 
axial-vector states can be categorized into the conventional axial-vector meson family. In this paper we carry 
out an analysis based on the mass spectra and two-body Okubo-Zweig-Iizuka-allowed decays. Besides testing 
the possible axial-vector meson assignments, we also predict abundant information for their decays and the 
properties of some missing axial-vector mesons, which are valuable for further experimental exploration of the 
observed and predicted axial-vector mesons. 

PACS numbers: 14.40.Be, 12.38.Lg, 13.25.Jx 


I. INTRODUCTION 

Among the light unflavored mesons listed in the Particle 
Data Group (PDG) [1], there are abundant light axial-vector 
mesons with a spin-parity quantum number J^ - 1^, which 
form a P-wave meson family. Usually, we adopt hi, bi, fi, 
and a I to express the corresponding states with the quantum 
numbers = O^fU-), U(l+^), 0+(U+), and r(U+), 

respectively. In Table I, we collect the experimental informa¬ 
tion on the observed hi, bi, fi, and ai states, as well as the 
corresponding resonance parameters and the observed decay 
channels. 

Facing so many axial-vector states in the PDG, we need to 
examine whether all of these states can be categorized into the 
axial-vector meson family, which is crucial for revealing their 
underlying structures. We also notice that most axial-vector 
states are either omitted by the PDG or are recent findings 
needing confirmation. Due to the unclear experimental status 
of light axial-vector states, we need to carry out a quantita¬ 
tive investigation of them, which would be helpful for further 
experimental studies, especially of those axial-vector states 
either omitted by the PDG or unconfirmed by other experi¬ 
ments. 

In this work, we carry out a systematic study of the axial- 
vector states by analyzing mass spectra and Okubo-Zweig- 
lizuka (OZI)-allowed two-body strong decay behaviors. Our 
investigations are based on the assumption that all of the ax¬ 
ial mesons can be explained within the conventional qq pic¬ 
ture. Comparing our numerical results with the experimental 
data, we can further test the possible assignments of the states 
in the axial-vector meson family. In addition, information 
on the predicted decays of the axial-vector states observed or 
still missing in experiments is valuable to further experimental 
study of axial-vector meson. 
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This paper is organized as follows. In Sec. 11, we present 
the phenomenological analysis by combining our theoretical 
results with the corresponding experimental data; the Regge 
trajectory analysis is adopted to study mass spectra of the 
axial-vector meson family and the quark-pair creation (QPC) 
model is applied to calculate their OZI-allowed strong decay 
behavior. Finally, the discussion and conclusion are given in 
Sec. III. 


II. PHENOMENOLOGICAL STUDY OF OBSERVED 
AXIAL-VECTOR STATES 

A Regge trajectory analysis is an effective approach to 
study a meson spectrum [31], especially a light-meson spec¬ 
trum. The masses and radial quantum numbers of light 
mesons with the same quantum number satisfy the following 
relation: 

+ ( 1 ) 

where Mo and M are the masses of the ground state and the 
corresponding radial excitation with radial quantum number 
n, respectively, denotes the slope of the trajectory with a 
universal = 1.25 + 0.15 GeV^ [31]. 

In Fig. 1 , we present the Regge trajectory analysis, in which 
we consider all of the axial-vector states listed in the PDG as 
shown in Table I. Besides the observed ones, we also predict 
some missing states and show them in Fig. 1. Additionally, 
we notice that there are two possible candidates for the ai 
meson with quantum number - 3^Pi, i.e., ai(1930) 

and ai(2095). On the other hand, both /i(1420) and /i(1510) 
can be an ss partner of /i(1285) by analyzing only the Regge 
trajectory. Thus, a further study of their strong decay behav¬ 
iors would help to test these possible assignments to the ob¬ 
served axial-vector states and could provide more predictions 
of the observed and still-missing axial-vector mesons, which 
are valuable for the future experimental exploration of axial- 
vector mesons. 

To obtain the decay behaviors of the axial-vector mesons, 
we adopt the QPC model, which was first proposed by Micu 
[32] and further developed by the Orsay group [33-37]. This 
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TABLE I: Resonance parameters and strong decay channels of the axial-vector states collected in the PDG [1]. The masses and widths are 
average values taken from the PDG. The states omitted from the PDG summary table are marked by a superscript tl, while the states listed as 
further states in the PDG are marked by a superscript b. 


State 

Mass (MeV) 

Width (MeV) 

The observed decay channels 

a, (1260) 

fli(1640)^ 
r(l++) 0,(1930)'' 

0,(2095)'’ 

0,(2270)'’ 

1230 ± 40 

1647 ±22 

1930/30 

2096 ± 17 ± 121 

2270+35 

250 ~ 600 

254 ± 27 

155 ±45 

451 ±41 ±81 

305+™ 

-’+'-’-40 

3,71 [2], Tip [3], 0-11 [4] 

37r [5], Tip [4, 6] , crn [5], /2(1270);r [5] 

3;r“ [7] 

idiCn^ [8] 

37rO [7] 

h,(1235) 

!+(!+-) h,(1960)'’ 

Z7,(2240)'’ 

1229.5 ± 3.2 

1960 ± 35 

2240 ± 35 

142 ±9 

345 ± 75 

320 ± 85 

IDTl [9-11] 

w;r0 [12] 

w;r0[12] 

,/'i(1285) 

/,(1420) 

0+(l++) /i(1510)'' 
/,(1970)^ 
/,(2310)^ 

1282.1 ±0.6 

1426.4 ± 0.9 

1518 ±5 

1971 ± 15 

2310 ±60 

24.2 ± 1.1 

54.9 ± 2.6 

73 ±25 

240 ± 45 

255 ± 70 

pV [13], qi^Ti [14-16], floTT [15-17], KKn [15, 16, 18] 

KKn [19, 20], KK'‘{%92)+c.c [18-20] 

^:^*(892)±c.c [21, 22], ;r+;r-p'[23] 

[24] 

[24] 

h,(1170) 

/ii(1380)^ 

0-(l+-) /ii(1595)^ 

;i,(1965)'’ 

h,(2215)'’ 

1170 ±20 

1386 ± 19 

1594 ± 15/^0 

1965 ±45 

2215 ±40 

360 ± 40 

91 ±30 

384±60|™ 

345 ± 75 

325 ± 55 

Tip [25-27] 

A'k'*(892)±c.c [21, 28] 

wq [29] 

U)q [30] 

UJq [30] 


model was widely applied to study the OZI-allowed two-body 
strong decay of hadrons [38-59]. In the following, we briefly 
introduce the QPC model. 

For a decay process A ^ B + C, we can write 


can be treated separately. In addition, i and j denote the color 
indices of a qq pair. 

By the Jacob-Wick formula [60], the decay amplitude is 
expressed as 


{BC\T\A) = H- , (2) 

where Pb(c) is a three-momentum of a meson B{C) in the rest 
frame of a meson A. A subscript Mj. (i = A, B, C) denotes 
an orbital magnetic momentum. The transition operator T is 
introduced to describe a quark-antiquark pair creation from 
vacuum, which has the quantum number 7^*' = 0 ^^, i.e., T 
can be expressed as 

T - -3y^(lm;l - m|00) j c/p3tfp45^(p3-i-P 4 ) 

{col% ^L(P3)<(P4), (3) 




^|A7T{2L+\) 
+ 1 


2 (L0;7MyJ74My,) 




and the general decay width reads 


r - 

4 ^ 


A J,L 


(5) 


where is the mass of an initial state A. We use the simple 
harmonic oscillator wave function to describe the space wave 
function of mesons, which has the following expression: 


which is constructed in a completely phenomenological way 
to reflect the creation of a quark-antiquark pair from vacuum, 
where the quark and antiquark are denoted by indices 3 and 
4, respectively. As a dimensionless parameter, y depicts the 
strength of the creation of qq from vacuum, where y = 8.7 
and 8.7/ VS [51] corresponds to the uuidd and ss creations, 
respectively. J/fm(p) = Ipl^frm(p) is the solid harmonic, (p, 
and u) denote the spin, flavor, and color wave functions, which 


'i'nImiR, p) = P)J//m(P), (6) 

where the concrete values of the parameter R involved in our 
calculation are given in Ref. [61] for the ground states. How¬ 
ever, its value is to be fixed for each excited state. 

With the above preparation, we further discuss the OZI- 
allowed decay behaviors of the axial-vector mesons. The al¬ 
lowed decay modes are listed in Tables II and III. 
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FIG. 1: (Color online) Regge trajectory analysis for ai, b\, fi, and hi with typical = 1.30 GeV^, 1.19 GeV^, 1.10 GeV^, and 1.19 GeV-, 
respectively, which can he covered by /r- = 1.25 ± 0.15 GeV- given in Ref. [31]. The experimental errors of discussed axial-vector states are 
given, which are taken from the PDG [1]. Here, o and • denote theoretical and experimental values, respectively. 


A. fli states 


The Regge trajectory analysis indicates that ai(1260) can 
be regarded as a ground state. The obtained total and partial 
decay widths of fli(1260) are listed in Fig. 2, which shows 
that np is the dominant channel. In Fig. 2, we give the partial 
decay widths of ai(1260) —» up from the S -wave and D-wave 
contributions. Flere, the D-wave/5 - wave amplitude ratio in 
the decay ai(1260) —> np is -0.248 with a typical value of 
R - 3.846 GeV * [61] in our calculation, which is compara¬ 
ble with the B852 data (-0.14 + 0.04 + 0.07) [4]. Our result 
also shows that fli(1260) ^ f^n is a subordinate decay mode 
with the partial decay width 1.82 MeV, which explains why 
there has been no evidence of ai(1260) ^ fon in experiments 
[3]. As shown in Fig. 2, the calculated total width can repro¬ 
duce the CMD2 data given in Ref. [62]. In addition, we also 
give some typical ratios relevant to the partial decay and to¬ 
tal widths together with the corresponding experimental data 
in Table IV. In summary, our results are comparable with the 
experimental values and support ai(1260) as a ground state in 
the a\ meson family. 

If fli(1640) is the first radial excitation of ai(1260), its 
decay behavior depending on the R value is shown in Fig. 
3. We use the experimental total width [5] and the ratio 
F(/2(1270);r)/F(cr7r) = 0.24+0.07 [5] to get R = (4.30 ~ 4.64) 



FIG. 2: (Color online) Total and partial decay widths of fli(1260) 
depending on R. Here, the dot-dashed line with the band is taken 
from the experimental data from Ref. [62]. The 5 -wave and O-wave 
contributions to the decay width of ai(1260) —> np are also given 
separately. All results are in units of MeV. 
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TABLE II; OZI-allowed two-body decay channels for ai and hi states marked by Here, p, ai, and 77 denote p(770), £j( 782), and 7/(548), 
respectively. The axial-vector states predicted by the Regge trajectory analysis are marked by a superscript tj. 


Channel 

01 (1260)oi(1640)o 1 (1930) oi (2095) a 1 (2270) 

Channel 

/7i(117O)/ii(138O)/7i(1595)/75(1780)/ii(1965) /75(2120)/!i(2215) /75(2340) 

Tip 

V 

V 

V 

V 

V 

np 

V 

V 

V 

V 

V 

V 

V 

V 

an 

V 

V 

V 

V 

V 

KK* 


V 

V 

V 

V 

V 

V 

V 

^/o 

V 

V 

V 

V 

V 

rjoj 


V 

V 

V 

V 

V 

V 

V 

^/i(1420) 


V 

V 

V 

V 

(oa 



V 

V 

V 

V 

V 

V 

;rp(1450) 


V 

V 

V 

V 

KKi(mO) 




V 

V 

V 

V 

V 

piii 


V 

V 

V 

V 

£J7/'(958) 




V 

V 

V 

V 

V 

77ao(980) 


V 

V 

V 

V 





V 

V 

V 

V 

V 

KK- 


V 

V 

V 

V 

poo(980) 




V 

V 

V 

V 

V 

;rZ7i(1235) 


V 

V 

V 

V 

;rp(1450) 




V 

V 

V 

V 

V 

^/2(1270) 


V 

V 

V 

V 

lf/fi(1400) 





V 

V 

V 

V 

^/i(1285) 


V 

V 

V 

V 

7fA'*(1410) 





V 

V 

V 

V 

pao(980) 



V 

V 

V 

7fii:*(1430) 





V 

V 

V 

V 

KKiiimO) 



V 

V 

V 

KK;(IA30) 





V 

V 

V 

V 

7/01(1260) 



V 

V 

V 






V 

V 

V 

V 

;rp(1700) 



V 

V 

V 

7/01(1420) 





V 

V 

V 

V 

KKiiniQ) 



V 

V 

V 

o-/ji(1170) 





V 

V 

V 

V 

KK^dAW) 



V 

V 

V 

;rp(1700) 





V 

V 

V 

V 

KK*(\A30) 



V 

V 

V 

pfl2(1320) 






V 

V 

V 

KKlimO) 



V 

V 

V 

£0/2(1270) 






V 

V 

V 

K-K" 



V 

V 

V 

£707(1420) 






V 

V 

V 

7/02(1320) 



V 

V 

V 

£0/i(1285) 






V 

V 

V 

era i(l 260 ) 



V 

V 

V 

p;r(1300) 






V 

V 

V 

crfl2(1320) 



V 

V 

V 

pai(1260) 






V 

V 

V 

p;r(1300) 




V 

V 

K'-Kiiino) 







V 

V 

7/oo(1450) 




V 

V 

fohi(mo) 







V 

V 

ojZ7i(1235) 




V 

V 

KK’‘d6m) 







V 

V 

p/7i(1170) 




V 

V 

£o/i(1420) 







V 

V 

po] (1260) 




V 

V 

K'KidAOO) 








V 

/C All (1270) 





V 










po2(1320) 





V 










pto(1420) 





V 










poo(1450) 





V 










^:a:*(1680) 





V 










7/'(958)ao(980) 





V 











GeV“". The main decay modes of ai(1640) are np, 7rp(1450), 
7r/2(1270), 7r/i(1285), and pat. Additionally, we provide fur¬ 
ther information on the typical ratios of ai(1640) decays in 


* Using the experimental total width [5], we find that overlap exists be¬ 
tween our theoretical and experimental results when taking R = 4.26 ~ 
4.92 GeV^'. Then, we can further constrain the R values by the ratio 
r(/2(1270)ff)/r(crff) = 0.24 ± 0.07 [5], where the constrained R = (4.30 ~ 
4.64) GeV“*, which is adopted to present other typical ratios of up 1640). 


Table V. 

There are two possible candidates for the second radial ex¬ 
citation of ai(1260). In the following, we discuss the decay 
behaviors of ai(1930) and fli(2095) by combining the corre¬ 
sponding experimental data. In Figs. 4 and 5, we present the 
R dependence of the decay behaviors of these ai’s, respec¬ 
tively. That is, the obtained total width of ai(1930) can be fit¬ 
ted with the data from Ref. [7] when R - 4.58 ~ 4.92 GeV ', 
while that of ai(2095) can overlap with the experimental data 
[7] when R = (4.78 ~ 5.16) GeV '. Thus, it is difficult to 
distinguish which ai is more suitable as a candidate for the 
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TABLE III: OZI-allowed two-body decay channels for bi and /i states marked by y. Here, p, a), and 77 denote p(770), w(782), and 77(548), 
respectively. The axial-vector states predicted by the Regge trajectory analysis are marked by a superscript \\. 


Channel 

*1(1235) *5(1640) *i(1960) *i(2240) 

Channel 

/i(1285)/i(1420)/i(1510)/f(1640)/f(1800)/i(1970)/f(2110)/J’(2210)/i(2310) 

TTO) 

V 

V 

V 

V 

;rao(980) 

V 

V 

V 

V 

V 

V 

V 

V 

V 

7rao(980) 

V 

V 

V 

V 

0-77 

V 

V 

V 

V 

V 

V 

V 

V 

V 

mi(1260) 


V 

V 

V 

;rai(1260) 


V 

V 

V 

V 

V 

V 

V 

V 

;ra2(1320) 


V 

V 

V 

KK* 


V 

V 

V 

V 

V 

V 

V 

V 

;rm(1420) 


V 

V 

V 

^2(1320) 



V 

V 

V 

V 

V 

V 

V 

;rao(1450) 


V 

V 

V 

;rao(1450) 




V 

V 

V 

V 

V 

V 

Tjp 


V 

V 

V 

7/0 




V 

V 

V 

V 

V 

V 

PP 


V 

V 

V 

PP 




V 

V 

V 

V 

V 

V 

KK* 


V 

V 

V 

(Jj(x) 




V 

V 

V 

V 

V 

V 

crp 


V 

V 

V 

ar]' 




V 

V 

V 

V 

V 

V 

y(958)p 



V 

V 

KKx(\21iS) 





V 

V 

V 

V 

V 

Pfo 



V 

V 

K^K- 





V 

V 

V 

V 

V 

mflo(980) 



V 

V 

07*1(1170) 






V 

V 

V 

V 

KKi(mO) 



V 

V 

^(95&)fo 






V 

V 

V 

V 

KKiilAOO) 



V 

V 

77/2(1270) 






V 

V 

V 

V 

KK;{1410) 



V 

V 

7fA-i(1400) 






V 

V 

V 

V 

KK;(1430) 



V 

V 

KK’-ilAlO) 






V 

V 

V 

V 

KK^imO) 



V 

V 

KK*{1430) 






V 

V 

V 

V 

o-/7i(1235) 



V 

V 

/:/f2*(1430) 






V 

V 

V 

V 

/CTC 



V 

V 

7/i(1285) 






V 

V 

V 

V 

0702(1320) 




V 

<r/2(1270) 






V 

V 

V 

V 

m;r(1300) 




V 

o-/i(1285) 






V 

V 

V 

V 

^r*^:i(1270) 




V 

0-/1 (1420) 







V 

V 

V 

77p(1450) 




V 

p*i(1235) 







V 

V 

V 

^:a:*( 1680 ) 




V 

77/1 (1420) 







V 

V 

V 

flo(980)/7i(1170) 




V 

0707(1420) 








V 

V 

P/2(1270) 




V 

7C7fi(1270) 








V 

V 

p/i(1285) 




V 

/:/C(1680) 









V 

P/i(1420) 




V 

/o/2(1270) 









V 

pii(1235) 




V 

/o/i(1285) 









V 

wai(1260) 




V 

oo(980)ai(1260) 









V 






flo(980)7r(1300) 









V 






00(980)02(1320) 









V 






77'(958)/2(1270) 









V 






PP(1450) 









V 






77'(958)/i(1285) 









V 






AT* a:, (1400) 









V 


second radial excitation of ai(1260) by studying only the total 
decay widths. Besides, we can learn from the Regge trajec¬ 
tory analysis that there is only one state for the state, 
and it is doubtful that both ai(1930) and a 1 (2095) exist, as 
mentioned in Ref. [7]. However, there exist different behav¬ 
iors of the partial decay widths of these ai’s. The ai(1930) 
mainly decays into final states up, 7rp(1450), and nb\{\2'i5). 


while the 7r/i(1285) and crn modes also have sizable contribu¬ 
tions. The decays of ai(1930) into /rA'p(1430), KK^{\A3Q), 
and /r*(896)/r*(896) have tiny decay widths, which are not 
listed in Fig. 4. As for ai(2095), its dominant decay channels 
are nb\{\21>5), np, and 7rp(1450) and are shown in Fig. 5. The 
other decay channels—like pao(980), 7rp(1700), 7r/i(1285), 
tt/o, and cr;7r—also have considerable contributions to the to- 
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TABLE IV: Some typical ratios of decay widths of ai(1260). The 
T{np)s^D) represent the S (D)-wave decay width of ai(1260) —> np. 

Our work Experimental data 

mvhyrZIi 086 0.60 [3] 

r((7rp)o)/r,-„„, 5.3 X 10-2 (1.30 ± 0.60 ± 0.22) x IQ-^ [3] 
T^^/Trotai 8.2 X 10-2 (18.76 ± 4.29 ± 1.48) x 10-2 [3] 
0-09_0.06 ±0.05 [1] 




R (GeV^) 


FIG. 3: (Color online) R dependence of the decay behaviors of 
ai(1640). Here, the dot-dashed line with the hand is the experimental 
total width from Ref. [5]. All results are in units of MeV. 


tal decay width. In Table VI, we also list some typical ra¬ 
tios relevant to their decays. We still need to emphasize one 
point. At present, ai(1930) and ai(2095) are not well estab¬ 
lished in experiments. The authors of Ref. [7] indicated that 
(22(1950) and (2i(1930) are not securely identified in mass and 
width, though some such contributions are definitely required 
[7]. However, when considering the Regge trajectory analy¬ 
sis, one finds that the 3^Pi state in the a\ meson family has a 
mass around 2000 MeV. The two unconfirmed (2i(1930) and 
(2 i( 2095) states could be candidates for the if’Pi state in the ai 
meson family, since their masses are close to that of the 3^Pi 
state in the ai meson family. Thus, the experimental study of 
the partial decay widths of (2 i( 1930) and (2 i( 2095) will help to 
reduce the two possible candidates—(2 i( 1930) and (21 (2095)— 
of the second radial excitation of the (2i(1260) to one. In 


the following, the experimental confirmation of (3i(1930) and 
(2 i( 2095) will be crucial for identifying the candidate of the 
3^Pi state in the ay meson family. If (3i(1930) and (2i(2095) 
cannot be established in experiments, we suggest an experi¬ 
mental search for (2i(3^Pi); the results for (2i(3^T’i) predicted 
in this work would be helpful for such a search. 


TABLE VI: Typical ratios for (2i(1930) and (2i(2095). The R ranges 
are (4.58 ~ 4.92) GeV-‘ and (4.78 ~ 5.16) GeV-‘ for ai(1930) and 
(2i(2095), respectively. 


Ratio 

ai(1930) 

fli(2095) 

^np/^Total 

0.151 -0.162 

0.139-0.176 

Pjrbi{n35)/^Tolal 

0.092 - 0.160 0.206 - 0.2542 

PirpinOOl/^Tolal 

0.005 - 0.024 0.039 - 0.0529 

^ crn!^ Total 

0.088 - 0.097 

0.058 - 0.073 

^7tp{ 1 450) I^T otal 

0.339 - 0.347 

0.189 - 0.253 

rj-i,,(1235)/rn-p(1450) 

0.271 - 0.462 

0.348 - 1.813 

r,ai(1260)/rA:A-*(892) 

0.629 -0.719 

1.141 - 1.742 

Ppiff„f2{m0) 

0.705 - 0.850 

0.188 - 0.451 

^' 7700 ( 980 )/r';rp( 1700) 

0.317 - 0.809 

0.239 - 0.279 

rA:A:](i400)/r,a2(i320) 0.508 ~ 0.553 

1.693 - 4.846 

F KK-i ( 1400) /F paa (980) 

- 

0.145 - 0.184 

F,ao(1450)/F)i:A:»(l430) 

- 

0.206 - 0.838 





FIG. 4: (Color online) R dependence of the calculated partial and 
total decay widths of the ai(1930). Here, the dot-dashed line with 
band is the experimental total width from Ref. [7]. All results are in 
units of MeV. 


TABLE V: Typical ratios of the decay widths of (3i(1640) corre¬ 
sponding to the R range (4.30 ~ 4.64) GeV-*. 


Ratio 

Value 

Ratio 

Value 

^ np!^ Total 

0.216 - 0.227 

Fj-p(1450)/FrDfa( 

0.473 - 0.474 

Fj-i,](1235)/Frorij/ 

0.014 - 0.059 

PkK’ IPa-n 

0.166 - 0.221 

F;r/ 2 ( 1270 )/Fpw 

0.523 - 0.855 F,jy|(i420)/Fi/|(i285) 0.089 - 0.094 

P’rfa 

0.166 - 0.221 




In Fig. 6, we discuss the decay behavior of aflllO) 
as the third radial excitation of (2i(1260). We find that the 
main decay mode includes the decay channels 7Tbi(l235), np, 
7rp(1450), and 7rp(1700). In addition, 7^7^‘(1410), p/ 2 i(l 170), 
/r/r*(1680), TTcr, and crai{l26Q) have important contribu¬ 
tions to the total decay width. p(22(1320), 77'(958)(2o(980), 
and K*Ki{l27Q) are subordinate decay modes, which are not 
shown in Fig. 6. In Table VII, we also list the typical ratios of 
the decays of the ai(2210). 
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FIG. 5: (Color online) R dependence of the calculated partial and total decay widths of ai(2095). Here, the dot-dashed line with the band is 
the experimental total width from Ref. [7]. All results are in units of MeV. 




R (GeV^) 


FIG. 6: (Color online) R dependence of the calculated partial and total decay widths of ai(2270). Here, the dot-dashed line with the band is 
the experimental total width from Ref. [7]. All results are in units of MeV. 


B. hi states 

The Regge trajectory analysis indicates that 7>i(1235), 
7>i(1960), and 7>i(2240) are the ground state, second radial ex¬ 
citation, and third radial excitation in the bi meson family, 
respectively. In addition, we also predict a missing 7>i(1640) 
as the first radial excitation. In the following, we study their 
decays. 

As for 7>i(1235), there are two allowed decay channels: nto 
and 7rao(980). The result shown in Fig. 7 shows that the ob¬ 


tained total width overlaps with experimental data from Ref. 
[63]. Since b\ ljji occurs via S and D waves, we obtain 
the D-wawefS -wave amplitude ratio of the bi —> u)n process, 
which is 0.465 in our work; this is consistent with the Crystal 
Barrel data (0.45 + 0.04) [10]. On the other hand, the decay 
channel nf^ has a partial decay width that is less than 1 MeV. 


As a predicted b^ state, 7>i(1640) has the decay behavior 
listed in Fig. 8, where we take the same R range as that 
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TABLE VII: Calculated ratios of the decays of fli(2270). Here, all 
the results correspond to the R range (5.12 ~ 5.32) GeV“*. 

Ratio Value Ratio Value 


r^rp/Er otal 
^Total 
En-pj7oo /Ero/fl/ 

^crn /E Total 


0.164 
0.247 • 
0.052 . 
0.064 


Ej-p(i450)/ErD(a/ 0.134 • 

Epa[(1260)/Ett,i,i(1235) 0.789 ' 
Es:A:,(i4oo)/Epi,3|^ 0.352 • 

Ep;,[(u70)/Eka-j(14oo) 0.573 • 


0.184 E,7/,(1285)/E.. 0.313-0.435 

0.264 T^j(^gg2)/^t]ai(l26Q) 0.313 ~ 0.487 
' 0.056 E;r/j(i420)/Epa„(980) 0.404 ~ 0.469 
0.070 E,,„„(980 )/Ei/2(1270) 0.532 ~ 0.612 
0.157 E,„3,i320)/E./„ 0.099 - 0.131 

0.926 E,,a„(i450)/E(„i|(i235) 0.236 - 0.273 
0.446 E,3 jj(1260)/Ea:a:»( 1680) 0.256 - 0.297 
0.639EK;7.(i4io)/E£roj(i260) 0.633 - 0.638 




4 4.3 4.6 4.9 

R(GeV’) 


4 4.3 4.6 4.9 

R(GeV’) 


FIG. 8: (Color online) R dependence of the calculated partial and 
total decay widths of 7>i(1640). All results are in units of MeV. 


130 
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120 
115 
110 
105 
100 

3 3.4 3.8 4.2 

R(GeV^) 

FIG. 7: (Color online) R dependence of the calculated total decay 
width of iii(1235). Here, the dot-dashed line with the band is the 
experimental total width from Ref. [63]. The total decay width is in 
units of MeV. 



for fli(1640). ^ Its main decay channel is too( 980), while 
77 (32(1320), pp, 7707(1420), KK*, and ljji also have consider¬ 
able contributions to the total decay width. The total decay 
width is predicted to be 200 - 232 MeV. Table VIII shows 
some ratios that are relevant to the decays of l7i(1640), which 
is valuable for further experimental searches for this axial- 
vector state. 


^ Since bi(\(AQ) is a predicted state, we take the same R range as that of 
fll(1640) to predict the decay behavior of l)i(1640). This treatment is due 
to the fact that i2i(1640) is the isospin partner of ai(1640), which has a 
similar R range. 


TABLE VIII: Typical ratios for decays of iii(1640) corresponding to 
R = 420~490GeV_^_ 

Ratio Value Ratio Value 

E™„(980)/Er<„., 0.352 - 0.368 Tkk^IT^„ 0.324 - 0.347 
E,p/E,j(j(i420) 0.164 - 0.263 pp 0.565 - 0.681 


Assuming that l7i(1960) is the second radial excitation of 
/7 i( 1235), we present its total and partial decay widths in Fig. 
9. Our calculated total width can cover the experimental data 
given in Ref. [12]. Its main decay channels are 7rao(1450), ttw, 
77ao(980), and 7707(1420), while the partial decay widths of the 
decay modes 77i3i(1260), prj, and 7702(1320) are also consid¬ 
erable. We also obtain some ratios of partial decay widths of 
Z7 i( 1960),which are listed in Table IX. 



FIG. 9: (Color online) R dependence of the calculated partial and 
total decay widths of /7i(1960). Here, the dot-dashed line with the 
band is the experimental total width from Ref. [12]. Since the width 
of the KK^ mode is tiny, we do not list its contribution here. All 
results are in units of MeV. 


In Fig. 10, we show the decay behavior of l7i(2240) as 
the third radial excitation of l7i(1235). Additionally, its main 
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TABLE IX: Obtained ratios for decays of ii(1960). All results cor- 


respond to R = 

4.66- 5.16 GeV-' 

1 


Ratio 

Value 

Ratio 

Value 

rtrag{9S0)/^Total 0.186 - 0.235 

^pp/^Total 

0.028 - 0.031 

f';rtj>i420 otal 

0.088 - 0.107 

^cjn/^Tolal 

0.077 - 0.162 

Tpp / f;ra2 (1320) 

0.572 - 0.624 

^KK'mTlI^tip 

0.648 - 0.736 

fp/o /f ;rfli(1260) 

0.029 - 0.030 TicK^iuoQi/t^KKiimo) 

0.249 - 0.316 


decay modes are loti, ;7rfa;(1420), nao(980), and :Tao(1450). 
Of course, the decay modes pp, pb\{l21>5), 7ra2(1320), and 
7rai(1260) also have obvious contributions to the total decay 
width. For the convenience of further experimental studies 
of this state, we provide information on typical ratios of the 
partial width of Z7i(2240) in Table X. 


TABLE X: Calculated ratios for Z7i(2240) corresponding io R = 
5.20 - 5.54 GeV-‘. _ 

Ratio Value Ratio Value 


FffflQ(980) /Fr Ola! 
Fjr(jJi 42 o /Fr otal 
Fm2(1320)/From/ 
F,p(i450)/Fro((i/ 
Total 

Fp/o/Fpi,,(i235) 


0.097 - 0.128 
0.068 - 0.071 
0.075 - 0.102 
0.055 - 0.064 


F 4 -flQ( 1450) /Fr otal 
^ ujn 1^ Total 
F^-aj (1260) 

Fpp/Fj'o^fl/ 


0.131 - 0.232 
0.179 - 0.199 
0.057 - 0.066 
0.050 - 0.062 


0.042 - 0.057 Fpy|(i285)/F2ja[(i260) 0.678 - 0.819 
0.112 - 0.173Fp/2(i270)/F*:s:[(i400) 0.254 - 0.317 


Fa:ki(1270)/Fk*( 1430) 0.233 - 0.383 Fp,,/(95g)/r7(2i:.(i4io) 0.364 - 0.381 
Fa:s;*(1410)/Fa-«:*( 892) 0.903 - 0.950 F^a„(980)/F(ja|(i260) 0.158 - 0.204 


C. f\ states 

When discussing /i states, we need to consider the admix¬ 
tures of the flavor wave functions \nh) - {\uu) + \dd))l V 2 and 
|iS>. /i(1285) and/i(1420)//i(1510) satisfy 


I/i(1285)) 


" cos (p - sin ' 

' \nn) ' 

J/i(1420)//i(1510))^ 


sin(;i cos(p 

, I'**) . 


where both /i(1420) and /i(1510) are partners of /i(1285). 
(We present their decay behaviors below.) (ft denotes a mixing 
angle. This mixing angle was determined in a phenomeno¬ 
logical way [64] and is given hy (j) — (20 - 30)°, which is 
consistent with (p - (24];^ 7 )° reported by the LHCb Collabo¬ 
ration [65] and cp - (21+5)° from the updated lattice QCD 
analysis [ 66 ]. When calculating the decays of /i(1285) and 
/i(1420)//i(1510), we take the LHCb value (p = 24°. 

In Fig. 1, we have predicted that /i(1640) is the first ra¬ 
dial excitation of /i(1285) and that /i(1800) is a partner of 
/i(1640); these two predicted axial-vector mesons are related 
by 


' l/i(1640))' 


" cos-sin 01 

" \nn) ' 

, l/i(1800)), 


sin 01 cos 01 



In addition, there exist relations among /i(1970), the pre¬ 
dicted /i(2110) and /i(2210), and /i(2310), i.e.. 


' l/i(1970))' 


" cos 02 - sin 02 ' 

" |nn) ' 

J/i(2110))^ 


sin 02 cos 02 , 



and 


' l/i( 2210 ))' 


" cos 03 - sin 03 ' 

" \nn) ' 

J/i(2310))^ 


sin 03 cos 03 

, I'**) . 


Here, the mixing angles cpj (i — 1,2,3) cannot be constrained 
by our analysis. In the following discussions, we take a typical 
value (pi - (p - 24° to give the quantitative results. 

As for /i(1285), we show its partial and total decay widths 
in Fig. 11, where the calculated total decay width is in agree¬ 
ment with the experimental data from Ref. [67]. However, 
we notice that the calculated branching ratio for total — 

0.67 - 0.68, corresponding io R - (3.00 - 4.00) GeV ', 
which is a little bit larger than (36 + 7)% listed in the PDG 
[1]. The PDG data also shows that the branching ratio of 
its decay rjjui can reach up to [ 1 ], which is the 

main contribution to the total decay width of /i(1285). In this 
work, we study the processes /i(1285) ^ t/ct —» rjiui and 
/i(I285) ^ 7rao(980) —» tjtttt, which can be calculated using 
the QPC model. Thus, the decay width of /i(1285) rjcr ^ 
TJTTTT can be written as [44] 


r(/i rj + cr ^ rj + tttt) 
1 


'^'■7 -272—;—FTT’ 

^ J4mi (r - + (m^Ta-y 




( 11 ) 


where the interaction of cr with two pions can be described by 
the effective Lagrangian 


-Ca 


■ ga-Cr{2TT^TT + TtV). 


( 12 ) 


The coupling constant go- - 2.12 - 2.81 GeV is determined 
by the total width Fo- = 400 - 700 MeV [1], and the decay 
width reads as 


_ glA^ [(r-( 2 m.)V]F 2 
8TTr 2 ^/r 


(13) 


where A - V2 and 1 for tt^tt^ and tt^tt^, respectively. 

The process /i(1285) —» too(980) —» tjtttt is calculated in a 
similar way, and the equation is given by 


r(/i ao + TT ^ TJ + tttt) 

^ (r - + {majattf ’ 


(14) 


where the decay width for (3o(980) ^ t/jt is 


rflo(980)-^;77r(^) 

_ gig [(>• - {m,, + Tnyy)(r - (m^ - 
8TTr 2 Vr 


( 15 ) 
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FIG. 10: (Color online) R dependence of the calculated partial and total decay widths of hi (2240). Here, the dot-dashed line with the band is 
the experimental total width from Ref. [12]. We do not present the K"Ki(inO) contribution since this decay has a tiny width. All results are 
in units of MeV. 
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FIG. 11: (Color online) R dependence of the total and partial decay 
widths of /i(1285). We also present the decay width of /i(1285) —> 
T/nn via the intermediate channels rjcr and naoi9SO) (green band), 
and only from the intermediate channel rjcr (pink band). Here, the 
experimental total width from Ref. [67] is denoted by the dot-dashed 
line with the band. All results are in units of MeV. 


where the coupling constant gag — 1.262 ~ 2.524 GeV 
is determined by the total width of ao(980) (ra„( 980 ) = 
50 ~ 100 MeV). The final result of the width of 

/i(1285) —> 7rao(980) —> rinn includes the contributions from 
both and rjn^TT^. 

The decay width of /i(1285) ^ rjTTTT via both the interme¬ 
diate rjcr and rao(980) channels and only the intermediate rjcr 
channel are shown in Fig. 11. In addition, the decay width 
of /i(1285) —> rinn from the intermediate 7rao(980) chan¬ 
nel is comparable with the corresponding experimental data 
[(16 + 7)%] in the PDG [1]. 

In the following, we discuss the decay behaviors of 


/i(1420) and /i(1510) as partners of /i(1285). As for 
/i(1420), the obtained total decay width can overlap with the 
DM2 result [68], as shown in Fig. 12. Its main decay channel 
is KK*. Thus, the present study of decay of /i(1420) sup¬ 
ports the prediction that /i(1420) is a partner of /i(1285). As 
for /i(1510), its partial and total decay widths are listed in 
Fig .13, which shows that the calculated total decay width is 
larger than the experimental data [68]. Thus, the possibility 
that /i( 1510) is a partner of /i(1285) can be excluded. 



R (GeV^) R (GeV^) 


FIG. 12: (Color online) R dependence of the total and partial decay 
widths of /i(1420). Here, the experimental total width from Ref. 
[68] is shown by the dot-dashed line with the band. All results are in 
units of MeV. 


In Figs. 14 and 15, we further illustrate the decay prop¬ 
erties of the two predicted states /i(1640) and /i(1800). In 
addition, we also list some of their typical ratios, which are 
weakly dependent on the R value (see Table XI), where we 
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R (GeV"^) R (GeV"^) 



FIG. 13: (Color online) Total and partial decay widths of /i(1510) 
as a partner of /|(1285) (first row) and as the first radial excitation of 
/i(1285) (second row). The experimental total width from Ref. [69] 
is denoted hy the dot-dashed line with the band. All results are in 
units of Me V. 


TABLE XI; Some obtained ratios relevant to decays of /i(1640) and 
/i(1800). All values correspond to the R range (3.60 ~ 4.40) GeV^'. 


/id 640) 

./'ll 1800) 

Fm,(i260)/FrD(a/ 0.400 ~ 0.440 
Fm2(i320)/FrD(a/ 0.114 ~ 0.312 
0.244 ~ 0.254 

TKK-ITTo,ai 0.026 - 0.284 

TppITro,.! 0.102- 0.290 
F(j£j/F,rni(i260) 0.254 - 0.665 
F/rA:i(i270)/FM* 0.088 - 0.141 
F™,(i320)/FrD(a/ 0.043 - 0.162 


take R — (3.60 ~ 4.40) GeV '. From Figs. 14 and 15 and Ta¬ 
ble XI, we can obtain information on the main decay modes 
and the resonance parameters of the two predicted /i mesons. 

As for/i(1510), there also exists another possible assign¬ 
ment, i.e., /i(1510) can be a radial excitation of/i( 1285) since 
the mass of /i( 1510) is close to that of the predicted/i( 1640). 
Here, we use the mixing angle expression 

|/i(1510)) = cos(pi\nh) - sin(;ii|ss), (16) 

which is the same as /i(1640). Thus, we also further illus¬ 
trate the decay behavior of /i (1510) as a radial excitation of 
/i(1285) (see Fig. 13). Under this assignment, the obtained 
total decay width can be fitted with the LASS data [69]. The 
KK* mode also has a large contribution to the total decay 
width. These facts indicate the possibility that /i(1510) is a 
radial excitation of /i(1285). 

In Fig. 16, we show the R dependence of the decay be¬ 
havior of /i(1970) as the second radial excitation of /i(1285). 
Its main decay channels are XX*(1410), 7rao(980), 7rai(1260), 
and KK*. As a partner of /i(1970), the predicted /i(2110) 
mainly decays into XXi(1270), XX*(1410), and KK* and has 
a large total decay width, as shown in Fig. 17. 

The third radial excitation of /i(1285) is still missing in ex¬ 
periments. In this work, we predict/i (2210), and we calculate 


FIG. 14: (Color online) R dependence of the total and partial decay 
widths of /i(1640). All results are in units of MeV. 



R (GeV"^) R (GeV"’) 

FIG. 15: (Color online) R dependence of the total and partial decay 
widths of /i(1800). All results are in units of MeV. 


its total and partial decay widths (see Fig. 18). As a part¬ 
ner of this predicted /i(2210), /i(2310) has the decay prop¬ 
erties listed in Fig. 19, in which the experimental width [24] 
is depicted by our calculation when taking R - (4.58 ~ 5.10) 
GeV '. Its main decay channels are XXi(1270), KK*(1680), 
XX*(1410), and XX*. 


D. hi states 

Similar to the fi mesons, the following study of hi states is 
relevant to the admixtures of the flavor wave functions nh and 
ss. As the ground states in the hi meson family, /ri(l 170) and 
























































12 


400 

300 

200 

100 

0 

50 
40 
30 
20 ' 
10 


Total Wi(«h 

20 

16 

12 

8 

4 

n 

^^/^o(1450) 

><r/7ta2(1320)/^ 

PR, 

1.5 

1.2 

0.9 

0.6 

0.3 

0 

Tiyi270) 

. oti(958kX 

SPEC 

: 

MT1260) 

2.5 

t^285) 

0.2 

, 


2 

^^^1270) 

0.16 

KK;(1430)y' 



/ 

0.12 


oyi2M) 

1.5 


0.08 

/ 

KK^(1270)^^^ 


/ 




0.5 


0.04 

KK.,(1400) 


0 

vy 

0 



4.7 4.9 5.1 5.3 5.5 4.7 4.9 5.1 5.3 5.5 4 

R (GeV"’) R(GeV-^) 


7 4.9 5.1 5.3 5.5 
R(GeV-^) 


contribution to the total decay width, which explains why the 
0)7] mode was found in Ref. [73]. As the predicted partner of 
/ii(1595), /ri(1780) dominantly decays into KK*, as presented 
in Fig. 21. 

Figure 22 presents the results for /ri(1965), where the cal¬ 
culated total decay width can overlap with the Crystal Barrel 
data [30] when R - (5.02 ~ 5.28) GeV^'. Its main decay 
channels are np, 7rp(1450), and 7rp(1700), while crh\{\\10) 
also provides a considerable value. As a partner of /2i(1965), 
/ 2 i( 2120 ) is predicted in this work, where its main decay 
modes are KK*, KK*{\A\Q), and /r/rg(1430) (see Fig. 23 for 
more details on its decay properties). 

The total and partial decay widths of h\{22\5) and its part¬ 
ner /2 i( 2340) predicted in this work are listed in Figs. 24 
and 25, respectively. The main decay modes of /2i(2215) and 
/2 i( 2340) can be found in Figs 24 and 25. 


FIG. 16: (Color online) R dependence of the total and partial decay 
widths of /i(1970). The experimental total width from Ref. [24] is 
denoted by the dot-dashed line with the band. All results are in units 
ofMeV. 


/2i( 1380) satisfy 


' |/ii(1170)> ' 


sin 01 cos 01 

" |n«) ' 

, -|/ii(1380)) ^ 


^ - cos 01 sin 01 ^ 

. I'Ss) . 


(17) 

where the mixing angle 0\ is introduced, the first line of this 
equation is adopted in this paper, and the second line is used 
in Ref. [70]. In Ref. [70], Cheng obtained 6 »i ~ 82.7°. The 
lattice QCD calculation indicates 6i - 86 . 8 ° [71]. In addition, 
6i - 85.6° was obtained in Ref. [72]. In our calculation, we 
present our result as 6\ - 85.6°. 

The obtained partial and total decay widths of /2i(1170) 
and /2 i( 1380) are shown in Fig. 20. Our results indicate that 
/2 i( 1170) and /2i(1380) are suitable candidates for the ground 
states in the hi meson family. Our result that /2i(1380) mainly 
decays into KK* is consistent with the experimental fact that 
/2 i( 1380) has a dominant ss component [21, 28]. 

According to the Regge trajectory analysis in Fig. 1, 
/2 i(1595), /2i( 1965), and /2i(2215) are the first, second and 
third radial excitations of / 2 i(l 170). Here, /2i(1595), /2i(1965), 
and /2 i( 2215) have the same flavor wave functions as / 2 i(l 170) 
in Eq. (17). The mixing angle 6i in Eq. (17) is replaced by 
6 * 2 , 03 , and 04 for the corresponding hi states. As for these 
higher radial excitations, the mixing angles 0, {i — 2,3,4) 
were not well determined. Thus, we take a typical mixing 
angle 0,- = 85.6° to discuss the decay behaviors of /2i(1595), 
/2 i( 1965), and/2i(2215). 

As for /2 1 ( 1595), we find that the obtained total decay width 
is much smaller than 384+60];™g MeV measured by the BNL- 
E852 Collaboration [29]. Thus, we suggest performing a pre¬ 
cise measurement of the resonance parameters of /2i(1595), 
which would be helpful in clarifying this discrepancy. The 
result shown in Eig. 21 indicates that np is a dominant decay 
mode of /2i(1595). In addition, /2i(1595) ^ ojtj has a sizable 


III. DISCUSSION AND CONCLUSION 

Although there are abundant axial-vector states in the PDG 
[ 1 ], the properties of the observed axial-vector states are still 
unclear. The present unsatisfactory research status of the ob¬ 
served axial-vector states stimulates us to systematically study 
them, which will be helpful for revealing their underlying 
structures. As a crucial step, we have studied whether the 
observed axial-vector states can be categorized into the axial- 
vector meson family. 

In this work, we have discussed the observed axial-vector 
states by assigning them as conventional states in the axial- 
vector meson family, where both the analysis of the mass 
spectra and the calculation of their two-body OZI-allowed 
strong decays have been performed. 

In our calculation using the QPC model, we took different R 
ranges to produce the total width of the discussed axial-vector 
states. Eor the ai and bi states, for example, we listed the 
obtained R values for different states (see Table XII for more 
details). We found that the corresponding R values become 
larger with an increase in the radial quantum number, which 
is consistent with our understanding, i.e., the size of higher ra¬ 
dial excitations is larger than that of lower radial excitations. 
Thus, our calculation can reflect this phenomenon, which pro¬ 
vides a test of the reliability of our calculation. In addition, we 
also noticed that states with the same radial quantum number 
in the ai and bi families have similar R ranges, which reflects 
the fact that the ai state is the isospin partner of the corre¬ 
sponding bi state. 

When we discussed the decay behaviors of higher radial ex¬ 
citations in the fi and hi meson families, we fixed the corre¬ 
sponding mixing angles to match the numerical results, which 
is due to the absence of a theoretical study of these mixing 
angles, and these mixing angles cannot be determined by the 
present experimental data [1]. However, for the ground states 
of fi and hi, the situation is totally different, where the corre¬ 
sponding mixing angles are fixed by experimental data. Thus, 
in this work we adopted a very simple and crude approach, 
i.e., we took the same value of the mixing angle for the ground 
and the corresponding radial excitations. We expect more ex- 
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FIG. 17: (Color online) R dependence of the total and partial decay widths of /i(2110). All results are in units of MeV. 



FIG. 18: (Color online) R dependence of the total and partial decay widths of /i(2210). All results are in units of MeV. 


perimental data on radial excitations in the fi and hi meson 
families. Then we can carry out further theoretical studies by 
considering the effect of the mixing angles. 

In summary, this phenomenological analysis not only tested 
possible assignments of the axial-vector states, but also pre¬ 
dicted abundant information about their partial decays, which 
is valuable for further experimental studies of the observed 
states. In addition, we have also predicted some missing axial- 
vector mesons, and have roughly determined their mass values 
and decay behaviors. We have also suggested an experimen¬ 
tal search for the missing states; the BESIII and COMPASS 
experiments will be good platforms to carry out the study of 
light hadron spectra. 
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FIG. 19: (Color online) R dependence of the total and partial decay widths of /i(2310). The experimental total width from Ref. [24] is denoted 
by the dot-dashed line with the band. All results are in units of MeV. 



3.2 3.4 3.6 3.8 4 4.2 



FIG. 20: (Color online) R dependence of the total decay width of 
/ii(1170) (top) and /2i(1380) (bottom). Here, the dot-dashed lines 
with the band are from Refs. [25, 28]. All results are in units of 
MeV. 
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FIG. 21: (Color online) R dependence of the total and partial decay 
widths of /ri(1595) (top) and /ri(1780) (bottom). All results are in 
units of Me V. 



FIG. 22: (Color online) R dependence of the total and partial decay 
widths of /ji(1965). The experimental total width from Ref. [30] is 
denoted by the dot-dashed line with the band. All results are in units 
ofMeV. 


TABLE XII: The obtained R value for these discussed a i and bi states 


in this work. 

state 


R (GeV-‘) 

state 


R(GeV-‘) 

a,(1260) 

I'F, 

3.846 

ii(1235) 

l‘Fi 

3.704 

ai(1640) 

2^P, 

4.30 ~ 4.64 

ii(1640) 

2 ‘Fi 


ai(1930) 

3^Pi 

4.58 ~ 4.92 

i|(1960) 

3‘F| 

4.66-5.16 

ai(2095) 

i^Pi 

4.78 ~ 5.16 




a, (2270) 

A^Pi 

5.12-5.32 

1^1(2240) 

4'F| 

5.20 - 5.54 
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R (GeV"') R (GeV"^) 


FIG. 23: (Color online) R dependence of the total and partial decay 
widths of /ji(2120). All results are in units of MeV. 



FIG. 24: (Color online) R dependence of the total decay width of h[(22l5). The experimental total width from Ref. [30] is denoted by the 
dot-dashed line with the band. All results are in units of MeV. 
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